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We report a measurement of the flux-averaged neutral-current elastic difi-erential cross section for 
neutrinos scattering on mineral oil (CH2) as a function of four-momentum transferred squared. It is 
obtained by measuring the kinematics of recoiling nucleons with kinetic energy greater than 50 MeV 
^^ ■ which are readily detected in MiniBooNE. This differential cross-section distribution is fit with fixed 

}_( ' nucleon form factors apart from an axial mass, Ma, that provides a best fit for Ma ~ 1.39±0.11 GeV. 

Cu ' Additionally, single protons with kinetic energies above 350 MeV can be distinguished from neutrons 

and multiple nucleon events. Using this marker, the strange quark contribution to the neutral- 
current axial vector form factor at Q^ = 0, As, is found to be As — 0.08 ± 0.26. 
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I. INTRODUCTION. 

Neutrino-nucleon neutral-current elastic (NCE) scat- 
tering is a unique and fundamental probe of the nu- 
cleon. NCE scattering on a nuclear target such as car- 
bon may be viewed as scattering from the individual nu- 
cleons but may also include contributions from collec- 
tive nuclear effects. This process should be sensitive to 
nucleon isoscalar weak currents as opposed to charged- 
current quasi-elastic (CCQE) scattering which interacts 
only via isovector weak currents. Therefore, the NCE 
process can be used to search for strange quarks in the nu- 
cleon which may show themselves via the isoscalar weak 
current. In addition, the NCE process offers a comple- 
mentary channel to CCQE to investigate any substantial 
collective nuclear effects in a nucleus such as carbon. 

The MiniBooNE experiment [i|,[3 searches for neutrino 
oscillations at Fermilab and NCE events account for 18% 
of the total neutrino sample collected. A large fraction of 
these neutral-current events are readily observable in the 
MiniBooNE detector which uses pure mineral oil (CH2) 
as a detector medium. The detector is predominantly 
a Cherenkov detector, however, fluors presented in the 
mineral oil produce a small amount of scintillation light 
well below Cherenkov threshold for relativistic charged 
particles (such as electrons, muons, protons, etc.). The 
absence of prompt Cherenkov light allows for the identi- 
fication and measurement of the recoiling nucleons pro- 
duced in neutrino NCE scattering. 

The main result presented in this paper is a high- 
statistics measurement of the flux-averaged differential 
cross section as a function of Q^ for NCE scattering on 
CH2 in MiniBooNE for 0.1 GeV^ < Q^ < igs GeV^. It 
is presented as scattering from individual nucleons both 
bound (in carbon) and free (in hydrogen). However, it is 
acknowledged that nuclear effects may well be important 
for full understanding of this data. 

A preliminary measurement of the MiniBooNE NCE 
differential cross section had been performed early on us- 
ing 6.57 X 10"'^^ protons on the neutrino production tar- 
get (POT) [3J. The results presented in this paper are 
based on the entire data set corresponding in neutrino 
run mode, for a total of 6.46 x 10^° POT and a different 
event reconstruction. MiniBooNE has also collected anti- 
neutrino data which will be the presented in subsequent 
works. 

To characterize the NCE from carbon in MiniBooNE, 
the relativistic Fermi gas (RFC) model of Smith and Mo- 
niz is used as described in Section IIIDl Within this 
model, two fundamental parameters are employed: the 
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nucleon axial mass, M^, and the strange quark contribu- 
tion to the axial form factor. As. It is acknowledged that 
nuclear effects in carbon may cause this axial mass to be 
an "effective" value and not the same as that for scatter- 
ing from free nucleons. The formalism for NCE scattering 
on free nucleons, which is basic to any model-dependent 
approach to characterize the NCE data, is summarized 
in Appendix B. 

While the generally accepted value for Ma is 1.026 ± 
0.021 GeV [j|, recent experiments [^-Q measuring 
CCQE from nuclear targets have found it useful to em- 
ploy values that are 20-30% larger to fit the Q^ (squared 
four-momentum transfer) dependence of their observed 
yields. It may be that this increased value of Ma should 
be understood not as the Ma obtained for free nucleons 
but rather as a parameterization of neglected nuclear ef- 
fects [^, llfll- Regardless, an extraction of Ma from NCE 
scattering offers a complementary test to the Ma deter- 
mined from CCQE. 

The strange quark contribution to the nucleon spin 
at Q^ — 0, As, can be extracted via the NCE process 
within an RFC model. The NCE differential cross sec- 
tion at low Q^ is sensitive to As [ll| for both neutrinos 
and antineutrinos. The BNL E734 experiment measured 
these processes and reported As [13 - [l4j . However, ratio 
measurements offer the possibility for extracting As at 
low Q^ with reduced systematic errors. For example, a 
measurement of iy-p — > vp)/{v^n — » ^,p), has been pro- 
posed by the FINeSSE experiment [i^|. In this paper, a 
ratio of {up ^ vp)/{vN -^ vN) at Q^ > 0.7 GeV is used 
to extract As, where A'^ is a neutron or proton. 

In the following sections the MiniBooNE experiment is 
described, including a description of the neutrino beam- 
line, detector, neutrino flux prediction, and the cross- 
section model used to predict the rates of different neu- 
trino interactions in the detector. In Sections IIIII and HVl 
methods and techniques used in the NCE analysis are 
presented together with various results, including the 
NCE differential cross-section, the NCE/CCQE ratio and 
measurements of Ma and As. Section FVl contains a sum- 
mary of the paper. Also, Appendices A and B have 
additional information useful for interpreting the results 
which is not included in the main text. 



II. MINIBOONE EXPERIMENT 

A. Neutrino Beamline. 

The Booster Neutrino Beamline (BNB) at the Fermi 
National Accelerator Laboratory uses a beam of pro- 
tons with momentum 8.89 GeV/c to produce an intense 
and almost pure beam of z/^ with an average energy of 
about 800 MeV. Protons are extracted from the Fermi- 
lab Booster in 1.6 /iS pulses with ~ 4 x 10^^ protons in 
each beam pulse. They are delivered onto a beryllium 
target, where a secondary beam of mesons is produced in 
p-Be interactions. Mesons are passed through a magnetic 
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FIG. 1: Neutrino flux at the MiniBooNE detector for different 
types of neutrinos as a function of their energy as reported 
in 0] and [11]. 



horn, a device which focuses positively charged particles 
and defocuses negatively charged particles. Mesons de- 
cay in an air-filled decay pipe producing a beam of neu- 
trinos. Using the magnetic horn increases the neutrino 
flux at the MiniBooNE detector by a factor of ~ 6. The 
details on the BNB components can be found in Rcf . T] . 



C. Detector. 

The MiniBooNE detector is situated 541 m from the 
Be target, under 3 m of overburden in order to reduce 
cosmic backgrounds. It is a 12.2 m diameter spherical 
steel tank, filled with mineral oil. The tank is divided 
into two optically isolated regions, a signal region being 
an inner sphere of radius 5.75 m, and a veto region that 
is an outer shell with a thickness of 0.35 m. Photomul- 
tiplier tubes (PMTs) are used to detect photons emit- 
ted by the charged particles which are produced in the 
neutrino interactions. Charged particles may emit both 
Cherenkov and scintillation light. Information from all 
PMTs is used to identify and reconstruct the products of 
neutrino interactions. 

A total of 1520 8-inch PMTs 0| are instrumented 
in the detector. There are 1280 PMTs attached to the 
spherical barrier from the inside, in the signal region, fac- 
ing toward the center of the tank and distributed approxi- 
mately uniformly. The remaining 240 PMTs are placed in 
the veto region, and are used to tag charged particles en- 
tering or leaving the tank. The veto PMTs are mounted 
back-to-back, tangentially to the optical barrier, in order 
to have as much veto view as possible. Details on the 
MiniBooNE detector can be found elsewhere Q- 



D. Cross-Section Model. 



Neutrino Flux Prediction. 



The neutrino flux at the detector is calculated via a 
GEANT4-based [165 MC beam simulation. The simula- 
tion includes a full beam geometry, specified by shape, 
location and material composition of the BNB compo- 
nents. The MC generates protons upstream of the target 
and propagates them through the target, generating and 
propagating products of p-Be interactions through the 
rest of the simulated BNB. In the f^ flux at the Mini- 
BooNE detector, 96.7% of neutrinos are produced via 
7r+ — !■ /i+ + i>^ decay. The tt"*" production double differ- 
ential cross section used in the beam MC is based on a 
fit to an external measurement from the HARP experi- 
ment on the same target and with the same proton beam 
energy as in the BNB [17]. 

The neutrino flux prediction for different types of neu- 
trino species is shown in Fig. [T] Flux tables are available 
in Ref. [18]. The 7r+ production contribution to the neu- 
trino flux uncertainty is about 5% at the peak of the 
flux distribution, increasing significantly at low and high 
neutrino energies. Other contributions to the flux error 
include uncertainties on other mesons production cross 
sections, the number of POT, and the horn magnetic 
field El, M. 



Neutrino interactions within the detector are simulated 
with the NUANCE-v3 event generator [2l|, where the 
relativistic Fermi gas model of Smith and Moniz [22| is 
used to describe NCE scattering. Fermi momentum for 
carbon is taken to be 220 ± 20 MeV and binding energy 
34 ± 9 MeV. 

The contribution from strange quarks to the vector 
and axial vector form factors is taken to be zero. The 
error on As id taken to be 0.1. The value of Ma used 
in the MC is different for the quasi-elastic (both neutral 
and charged current) scattering on carbon and hydrogen: 
for scattering on carbon Ma — 1.230 ±0.077 GeV is used 
(as measured from the CCQE channel in MiniBooNE [5]), 
while for scattering on hydrogen Ma = 1.13±0.10 GeV is 
used (which is the average between the values measured 
by the deuterium-based scattering experiments and Mini- 
BooNE). 

For resonant pion production, the Rein and Sehgal 
model 23] is used. In the few GeV range, such pro- 
cesses are dominated by the A(1232) resonance, although 
contributions from higher mass resonances are also in- 
cluded in the MC. The values Ml"" = 1.10 ± 0.27 GcV 
for both charge current (CC) and NC single pion events, 
and A/^'^ = 1.30 ± 0.52 GeV for both CC and NC mul- 
tiple pion events are used. 

Intranuclear final state interactions (FSI) inside the 
carbon nucleus are modeled in NUANCE using a binary 
cascade model [21], where the scattered hadrons are prop- 
agated through the nucleus, which is simulated based on 



models of nuclear density and Fermi momentum. Due to 
FSI, a NCE interaction may produce more than one final 
state particle (other than the neutrino). For NCE scat- 
tering on carbon, the probability of producing multiple 
nucleons is ^ 26% integrated over the MiniBooNE flux, 
according to NUANCE. Also, a NC pion event might not 
contain any pions in the final state as the pion can be ab- 
sorbed in the carbon nucleus or the baryonic resonance 
re-interacts without decaying. These are the dominant 
mechanisms by which NC pion events can become back- 
grounds to this analysis. The probability that a pion is 
absorbed is ~ 20% in carbon for MiniBooNE energies, 
according to the NUANCE simulation. The intranuclear 
pion absorption cross section is assigned a 25% uncer- 
tainty based on external pion-carbon data 12414 26l |. and 
A7V — )■ NN interactions are assigned a f 00% uncertainty. 

Neutrino interactions outside the detector, in the sur- 
rounding dirt or in the detector material (referred to as 
"dirt" background henceforth), are simulated the same 
way as the in-tank interactions but with a cross section 
reweighted according to the density of the material rela- 
tive to that of the mineral oil. 

Particle propagation in the detector is modeled using 
a GEANT3-based [Hj MC with GCALOR [28,] hadronic 
interactions simulating the detector response to particles 
produced in the neutrino interactions. GCALOR was 
chosen over GFLUKA 29] simulation, as it p rovides a 
better model of 7r+ interactions on carbon [301 ■ 



III. NEUTRAL-CURRENT ELASTIC 
ANALYSIS. 

A. Event Reconstruction. 

Event reconstruction in MiniBooNE is based on find- 
ing a set of parameters (position, time, direction and 
energy ~ where applicable) which maximizes the event 
likelihood using the charge and time information from 
all PMTs. Each event is reconstructed using some com- 
bination of six different event hypotheses - single proton 
(NCE- like), single muon (i/^ CCQE-like), single electron 
(I'e CCQE-like), single tt" (NC tt" production-like), muon 
and 7r+ with the same vertex (i/^ CC tt"*" production- 
like), and muon and tt*' with the same vertex (i/^ CC 
tt" production-like). A charge-time likelihood minimiza- 
tion method |3l| is used to obtain the best estimate of 
the kinematic observables in each event hypothesis. Un- 
der the NCE hypothesis, each event is assumed to be a 
point-like proton with Cherenkov and scintillation light 
emission profiles determined from the MC. The output 
variables from these event reconstructions (such as likeli- 
hood ratios between two different event hypotheses) allow 
for particle identification. 

The resulting position resolution is ~ 0.75 m for proton 
events in the detector and ~ 1.35 m for neutrons, with 
an energy resolution of ~ 20% for protons and ^ 30% for 
neutrons. For protons above Cherenkov threshold, the 
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FIG. 2: Fraction of prompt hits versus the total number of 
tank PMT hits for beam unrelated data and NCE MC events 
reconstructed under an electron hypothesis. The error bars 
correspond to the RMS of the distributions. 



direction resolution is ~ 10°. 

Light emission properties of protons differ from those 
of other charged particles in the detector, allowing for 
their particle identification. For instance, protons differ 
from electrons in terms of the fraction of prompt light 
emitted, defined as the fraction of PMT hits with a cor- 
rected time between —5 ns and 5 ns, as illustrated in 
Fig. m where the corrected time is the time difference 
between the PMT hit time and the reconstructed event 
time, with light propagation time from the reconstructed 
vertex to the PMT also taken into account. 

Being neutral particles, neutrons themselves do not 
cause light emission in the detector. However we 
may detect them through their subsequent strong re- 
interactions, in which usually energetic protons are pro- 
duced. Because we detect NCE neutrons only through 
secondary protons, they are virtually indistinguishable 
from NCE proton events. 

For MiniBooNE NCE interactions, the total charge on 
all PMTs is proportional to the sum of kinetic energies 
of all final state nucleons that are produced in the inter- 
action, which is referred to throughout this paper as T. 
ft is important to understand that the nucleon kinetic 
energy measured this way is different from the one deter- 
mined from the track-based reconstruction used in the 
SciBooNE [32] and BNL E734 [13] experiments. In that 
case, the reconstructed proton track length is propor- 
tional to the kinetic energy of the most energetic proton 
produced in the event. Also the particle identification 
in MiniBooNE is based almost entirely on the properties 
of the measured Cherenkov ring (such as ring sharpness, 
charge and time likelihoods), whereas the track-based ex- 
periments mostly use the particle's energy loss along the 
track. 



B. Event Selection. 

The following set of selection criteria (cuts) are applied 
to the full MiniBooNE data set to select the NCE sample: 

1. Only 1 subevent to ensure the event is NC and 
includes no decaying particles (e.g., /i decay). A 
subevent is a cluster of at least 10 tank PMT hits 
for which there is no more than 10 ns between any 
two consecutive hits. For example, a CCQE event 
typically contains two subevents; the first subevent 
is associated with the outgoing muon, while the sec- 
ond is associated with the subsequent decay elec- 
tron [3. 

2. Number of veto PMT hits is less than 6 to remove 
events exiting or entering the detector, VHits < 6. 

3. Number of tank PMT hits is greater than 24 to 
ensure a reliable reconstruction, THits > 24. 

4. Beam time window cut in order to consider only 
events time-coincident with the neutrino beam. 

5. Reconstructed proton energy of T < 650 MeV 
(above which the signal to background ratio de- 
creases significantly). 

6. Log-likelihood ratio between electron and proton 
event hypotheses of \n{£e/^p) < 0.42. The pur- 
pose of this cut is to eliminate beam- unrelated elec- 
trons from cosmic-ray muon decays. Even though 
the variable shown in Fig. [2] provides a good sepa- 
ration between proton and electron events, it was 
found that the log-likelihood ratio is most effective 
in capturing all differences between the NCE sig- 
nal and this background. The \n{Ce/Cp) variable 
and the value of the cut are shown in Fig. |3] for 
simulated NCE and the beam unrelated data. The 
beam unrelated data events are mostly muon decay 
(Michel) electrons. 

7. A fiducial volume cut, defined as follows: 



^fiducial [-^ ) 



R<4.2m if T< 200 MeV, 
i?<5.0m if T> 200 MeV. 



A tighter fiducial volume is required at low energies 
to reduce the dirt background. 

A total of 94,531 events pass the NCE cuts resulting 
from 6.46 x lO^o POT. This is the largest NCE event 
sample collected to date. The efficiency of the cuts is es- 
timated to be 35%; a large portion of which stems from 
the fiducial volume cut. We consider all NCE events 
with original vertices inside the detector as signal. The 
predicted fraction of NCE events in the sample is 65%. 
The remaining 35% of events are backgrounds of different 
types: 15% are NCE-like backgrounds, 10% dirt events, 
and 10% other backgrounds (of which only 0.5% are beam 
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FIG. 3: Log-likelihood ratio between electron and proton 
event hypotheses for MC-generated NCE scattering events 
and beam unrelated data. Both histograms are normalized to 
unit area. Events with \n{Ce/Cp) < 0.42 are selected for the 
analysis. 



unrelated). The reconstructed nucleon kinetic energy 
spectrum for selected NCE events with a uniform fiducial 
volume cut is shown in Fig. 2] along with the predicted 
background contributions. 

The NCE-like background consists of NC pion produc- 
tion channels with no pion in the final state (i.e. the pion 
is absorbed in the initial target nucleus through FSI). In 
this case, the final state particles for these events are 
solely nuclcons. In MiniBooNE, this is indistinguishable 
from the final state produced in NCE events, and hence 
why these events are referred to as NCE-like background. 
The NCE-like background contributes mostly at interme- 
diate energies, 200 MeV < T < 500 MeV. 

The dirt background is an important contribution to 
the NCE data sample at low energies, most significantly 
below 200 MeV. This background is due to nucleons 
(mainly neutrons) which are produced in neutrino in- 
teractions outside of the detector, penetrating into the 
detector without firing enough veto PMTs. Dirt events 
are challenging to simulate in the MC because they oc- 
cur in various media that have not been studied in detail 
(in the soil, detector support structures, etc.). However, 
we directly constrain this background using MiniBooNE 
data. Dirt events can be isolated from in-tank interac- 
tions using their distinct kinematics: dirt events are pref- 
erentially reconstructed in the most upstream {Z < m) 
and outer regions of the detector with relatively low en- 
ergies (small values of T). The dirt energy spectrum is 
measured by fitting dirt-enriched samples in the variables 
Z, R and T, as explained in detail in Appendix A. 

Other backgrounds are mainly charged-current chan- 
nels, but also include neutral current pion production, 
beam unrelated, and anti-neutrino NCE events. These 
backgrounds become relevant at high reconstructed en- 
ergies, mainly above 400 MeV. 
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the interaction assuming the target nucleon is at rest. In 
this case we define: 
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FIG. 4: Reconstructed nucleon kinetic energy spectra for the 
data and MC after the NCE event selection and a uniform 
fiducial volume cut of _R < 4.2 m are applied. All MC dis- 
tributions are normalized to the number of protons on target 
(POT). 



C. Unfolding. 

After subtracting the beam unrelated and dirt back- 
grounds, the reconstructed energy spectrum for the data 
is multiplied by the signal fraction, which is the number 
of NCE events divided by the total number of neutrino in- 
duced in-tank events according to MC prediction, bin by 
bin. The obtained NCE reconstructed energy spectrum 
is then corrected for detector resolution and efficiency ef- 
fects that distort the original spectrum. The correction 
is applied by unfolding the distribution using a method 
based on Bayes' theorem 



The unfolding matrix applied to the reconstructed 
spectrum is calculated for NCE events in the MC us- 
ing a migration matrix of the true nucleon kinetic energy 
(the sum of the kinetic energies of all nucleons in the 
final state) versus the reconstructed nucleon kinetic en- 
ergy. This method gives a well-behaved but biased solu- 
tion, which depends on the original MC energy spectrum. 
The error due to the unfolding bias is estimated by it- 
erating the unfolding procedure for each MC variation, 
where the new MC energy spectrum is replaced by the 
unfolded energy spectrum. The details of the unfolding 
procedure used in the measurement and the error esti- 
mation associated with its bias can be found in Ref.(3l|. 
The detector resolution and cut efficiency effects can be 
seen by comparing the reconstructed energy spectrum be- 
fore unfolding (Fig. |4|) and the final cross-section result 
(Fig. [SJ, which we shall discuss later. 



D. Neutral— Current Elastic Flux- Averaged Cross 
Section. 

For each NCE event, Q^ can be determined by mea- 
suring the total kinetic energy of outgoing nucleons in 
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2mNT — 2mN N^ Ti, 



where ttijv is the nucleon mass and T is the sum of the 
kinetic energies of the final state nucleons. MiniBooNE 
reports the NCE differential cross section as a function 
of of this variable, Qqe- 

The MiniBooNE NCE differential cross-section is less 
sensitive to FSI effects than those measured by tracking 
detectors, such as BNL E734 12] . In case of a FSI, where 
no pions in the final state have been produced, the en- 
ergy transferred by the neutrino may be divided among 
several outgoing nucleons, but the total energy released 
in the MiniBooNE detector stays roughly the same due 
to energy conservation. Track-based detectors measure 
Q^ by the proton track length and its angle with respect 
to the beam direction, which are kinematic observables 
of the most energetic proton produced in the NCE event. 
In that case, FSI may have large effects on the kinemat- 
ics of individual outgoing nucleons, including the most 
energetic nucleon. Of course, there are still some FSI 
interactions producing final state pions which must be 
modelled and which can affect MiniBooNE NCE cross- 
section measurement. 

The resulting NCE flux-averaged differential cross sec- 
tion on CH2 is shown in Fig. [5] as a function of Qqe- The 
predicted distribution of the NCE-like background, which 
has been subtracted along with the rest of backgrounds, 
is also shown in the figure. The NCE scattering is a sum 
of three different processes: scattering on free protons in 
hydrogen, bound protons in carbon, and bound neutrons 
in carbon. A detailed description of the contributions 
of each of these processes to the total MiniBooNE NCE 
cross section is given in Appendix B. 

Systematic uncertainties and their contribution to the 
total error have been studied. The normalization error 
can be represented by a single number, shown in Ta- 
ble m The largest systematic error in the NC analysis, 
the optical model, arises from the uncertainty on both 
Cherenkov and scintillation light production by charged 
particles in mineral oil. Cherenkov light production has 
been calibrated on cosmic muons, Michel electrons and 
neutral pions [2|- An error of 20% on this amplitude of 
scintillation light generated by sub-Cherenkov particles 
in the detector has been assigned through a combination 
of benchtop measurements and in situ calibrations (see 
Appendix C of Ref ^T] for details) . 



E. Neutral— Current Elastic to Charged— Current 
Quasi-Elastic Cross-Section Ratio Measurement. 

Given that MiniBooNE measures a CCQE differen- 
tial cross section [l^ that is '^ 30% higher than naive 
expectations from the relativistic Fermi gas model t22|, 
it is interesting to compare those results with our NCE 
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FIG. 5: The MiniBooNE NCE {uN — >■ uN) flux-averaged differential cross section on CH2 as a function of Qqe = 2niN X^^ Ti, 
where we sum the true kinetic energies of all final state nucleons produced in the NCE interaction. The blue dotted line is the 
predicted spectrum of NCE-like background which has been subtracted out from the total NCE-like differential cross section. 



Type of error 


Error 


Value (%) 


Flux 


Flux 


6.7 


Cross section 


NC tt" production yield 
In-tank background cross section 
Dirt background 


0.5 
6.3 
1.0 


Detector 


Discriminator threshold 
Optical model 
Charge-time PMT response 
Hadronic interactions 


0.6 

15.4 

2.1 

0.5 


Total 




18.1 



TABLE I: Individual error contributions to the total inte- 
grated normalization uncertainty on the MiniBooNE mea- 
sured NCE cross section. The statistical error of 2.5% is not 
included in the total normalization error. 



measurement. To facilitate such a comparison and, at 
the same time, reduce flux uncertainties, we extract the 
NCE/CCQE ratio as a function of Qqe- In the case of 
CCQE, Qqe has been defined from the outgoing muon 
kinematics only, assuming a stationary neutron target 
(see Ref. pS| for details). It should be pointed out that 
a significant difference exists in how these cross sections 
are measured in MiniBooNE. As explained earlier, the 
NCE cross section is calculated from the measured total 
kinetic energy of final state nucleons and is mildly sen- 
sitive to FSI. Whereas the CCQE is calculated entirely 
from the reconstructed muon and is not sensitive to FSI. 



The measured ratio is shown in Fig. [5] together with 
the NUANCE MC prediction. The data/MC agreement 
is reasonable within errors. 



Adding the MC NCE-like background prediction to the 
numerator and the MC CCQE-like background predic- 
tion to the denominator produces a NCE-like to CCQE- 
like differential cross-section ratio, which is additionally 
shown in Figj?) This is an even more model independent 
measurement, where we do not have to rely on model- 
ing of both NCE-like and CCQE-like backgrounds and 
claim them as a part of the signal. 



The measured NCE/CCQE ratio is consistent with 
that predicted by the MC. This is an important point 
when considering possible explanations of the larger than 
predicted value of the CCQE cross section. The predicted 
MC ratio is chosen for two values of Ma and k: one_is 
with AIa = 1.23 GeV and n = 1.022 as measured in [5|, 
the second with Ma — 1.35 GeV and k = 1.007 is a from 
a more recent MiniBooNE CCQE result [l^, where n is 
a Pauli blocking scaling factor parameter. 



There is some disagreement between data and MC for 
the NCE-like/CCQE-like ratio above Q^^^ > 1-0 GeV^ 
but this is where the NCE-like backgrounds (predomi- 
nantly NC pion channels with pion absorption) become 
a significant fraction of the signal. 
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FIG. 6: MiniBooNE NCE/CCQE cross section ratio on CH2 
as a function of Qqe- The NUANCE MC prediction is plotted 
for two assumptions of parameters: black solid line for Ma = 
1.23 GeV and k. = 1.022_J1, and blue dotted line for Ma = 
1.35 GeV and k — 1.007 [ly|. As is assumed to be zero in both 
of these cases. Both cross sections in the ratio are per target 
nucleon - there are 14/6 times more target nucleons in the 
numerator than in the denominator. The error bars represent 
both statistical and all systematic uncertainties (excluding 
flux errors) taken in quadrature. 
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FIG. 7: MiniBooNE NCE-like/CCQE-like cross-section ratio 
on CII2 as a function of Qqe with total error. The color code 
for the histograms is the same as in the Figure [6l 



where gA = 1.2671 is measured precisely from neutron 
beta decay |34l |. 

To first order, the MiniBooNE i^N — > i^N cross section 
is not sensitive to As, as the linear term in As nearly 
cancels, while the quadratic term in As remains, but is 
small if I As I ^ gA- However, these data are still useful 
for probing Ma- 

A x^ goodness of fit test is performed to find the set of 
Ma and k parameters, described earlier in Section B.2, 
that best matches data. Varying values of Ma and k 
in the MC model results in different reconstructed NCE 
energy distributions. For each set of Ma and n values, a 
X^ is calculated using the full error matrix. The full er- 
ror matrix in this case also includes an additional uncer- 
tainty on As of 0.1. The reconstructed energy spectrum 
for data and MC for different values of Ma and k (as 
measured from MiniBooNE CCQE data, as well as the 
average values prior to MiniBooNE) are shown in Fig. [S] 
The reconstructed energy distribution of the NCE sam- 
ple has a negligible sensitivity to k; however, the higher 
values of Ma seem to better describe MiniBooNE NCE 
data. 

The NCE data can also be directly fit to indepen- 
dently extract information on Ma and As. Because of 
the quadratic term in As in the NCE differential cross- 
section expression, the shape of x^ slightly depends on 
the value of As. Assuming As = 0, the la allowed region 
of Ma from the MiniBooNE NCE sample yields: 



Ma^ 1.39 ±0.11 GeV, 

with x'Lin/DOF = 26.9/50. Using As = -0.2 (which 
roughly corresponds to the value obtained by the BNL 
E734 experiment [13) yields Ma = 1.35±0.11 GeV with 
xLin/DOF = 24.9/50. The results from the Ma fit to 
the NCE data using an absolute (POT) normalization 
agree well with the shape-only fit results from the Mini- 



BooNE CCQE data 



map ' 

m 



IV. MEASUREMENT OF As USING A HIGH 
ENERGY PROTON-ENRICHED SAMPLE. 

A. High Energy Proton-Enriched Sample. 



Axial Vector Mass Measurement Using the 
NCE Cross Section. 



In the Appendix B.l the expression for NCE differ- 
ential cross section on free nucleons is described. From 
there, one can see that the NCE cross section is sen- 
sitive to the axial form factor. In fact, at low Q^, 
da/dQ^ - (F|)2(g2), where 

FliQ') = i(5AT3 - As)/ (1 + Q'/MlY , 



As mentioned above, the As sensitivity of the NCE 
sample comes down to the possibility of distinguishing 
proton from neutron events. In particular, the ratio of 
vp -^ vp to vn — > vn is most sensitive to As in addition 
to reducing the systematic error. One should mention 
that even though FSI effects may be significant in the 
NCE cross section, they become negligible when using 
such ratios |35| 

In MiniBooNE, NCE scattering from a neutron can 
only be detected when that neutron has a further strong 
interaction, usually with a proton. At low energies, 
these cannot be distinguished from NCE scattering from 
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FIG. 8: x^-tests performed on the NCE reconstructed nucleon 
kinetic energy distribution for MC with different Ma values 
of 1.35, 1.23, and 1.02 GeV. The ^ values are 27.1, 29.2 
and 41.3 for 49 degrees of freedom (DOF), respectively. The 
distributions are normalized to POT. 



a proton. When scattering kinematics produces a pro- 
ton above Cherenkov threshold, it is distinctive since the 
secondary interaction required to detect a neutron rarely 
produces an above-Cherenkov-threshold proton. We thus 
select a special class of NCE protons, with only one pro- 
ton in the final state whose energy is above Cherenkov 
threshold. These single proton events will be used for a 
As measurement. 

It should be noted that events with multiple nucle- 
ons in the final state (both NCE neutrons and NCE-like 
background) have multiple protons produced in the reac- 
tion. At kinetic energies above 350 MeV, single proton 
events produce on average a higher Cherenkov light frac- 
tion than multiple proton events. In addition, these two 
classes of events differ in the kinematics of the outgoing 
nucleon, with single proton events being more forward- 
going. 



FIG. 9: The "Proton/Muon" cut based on the ln{Cf,/Cp) 
variable, which is the log-likelihood ratio of muon and proton 
event hypotheses. The left part of the cut is used for the 
analysis. The histograms are normahzed to POT. 
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FIG. 10: The "Proton/Neutron" cut based on the fraction of 
prompt hits with the corrected time between ns and 5 ns. 
The right part of the cut is used for the analysis. The his- 
tograms are normalized to POT. 



B. Ratio of i/p —5- jyp to i/N 



vN. 



In order to measure As, the ratio of vp -^ vp to 
vN -^ vN as a function of the reconstructed nucleon 
kinetic energy from 350 MeV to 800 MeV is used. Ad- 
ditionally, by taking the ratio, several sources of sys- 
tematic uncertainty are reduced. The denominator of 
this ratio are events with NCE selection cuts described 
in Section IIIIBl but with the energy cut (5) replaced 
with 350 MeV < T < 800 MeV, and an additional "Pro- 
ton/Muon" cut based on ln(£^/£p) (the log-likelihood 
ratio between muon and proton hypotheses) in order to 
reduce muon-like backgrounds that dominate the high 
visible energy region (Fig. |9l). After these cuts, the 
vN — > vN data sample includes 24,004 events with the 
following predicted channel fractions: 45% NCE, 26% 
NCE-like background, 3% dirt background, and 25% 



other backgrounds. 

For the numerator of the ratio, two more cuts are ap- 
plied in addition to the ones used for the denominator. 
The first of them, the "Proton/Neutron" cut, is a vari- 
able based on the fraction of prompt light (as described in 
Section III.B, but with a corrected time between ns and 
5 ns). This cut is to increase the single proton event frac- 
tion in the sample. The described variable distribution 
and the value of the cut are shown in Fig. [10] for both 
single proton and multiple proton events. As one can 
see, the cut reduces multiple proton events, which have 
less Cherenkov light than single proton events. Finally, 
we apply a cut on the angle between the reconstructed 
nucleon direction and the incident beam direction. As 
shown in Fig. Ill[ single proton events are mostly for- 
ward going; we thus require Op < 60°. The final vp — >■ vp 
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FIG. 11: The angle between the reconstructed nucleon direc- 
tion and the incident beam direction for data and MC after 
the NCE, "Proton/Muon" , and "Proton/Neutron" cuts. All 
distributions are POT-normalized. Events with 6p < 60° are 
used in the analysis. 
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FIG. 12: The ratio of up — >• vp/vN — >• vN as a function of 
the reconstructed energy for data and MC with As values as 
labeled. 



vp — !■ vp/vN — >■ vN ratio yields: 



As == 0.08 ± 0.26, 



(1) 



with xliin/DOF = 34.7/29. Using Ma = 1.23 GeV 
yields As = 0.00 ±0.30 with xfn,n/DOF = 34.5/29. The 
result is consistent with the BNL E734 measurement [12] . 

One needs to comment on the implications of Fig- 
ures [31 [Til and [121 on the As result. From Fig. [4l one 
can see that the MiniBooNE MC overpredicts the to- 
tal number of events passing the NCE selection cuts 
at high reconstructed energies {T > 250 MeV) by as 
much as 40%. The NCE proton-enriched sample was 
obtained for events with high reconstructed energies, 
350 < T < 800 MeV. For these events, looking at the 
6p distribution in Fig. [TTl it seems that the entire dis- 
agreement between data and MC at these energies comes 
from the forward-going events, vp — > vp. Clearly there 
is a deficit of vp — >■ vp in the data, which in principle 
implies that this might be due to positive values of As. 
However, Fig. [T^and the result in Eq.(|T]) shows that it 
is consistent with zero. This may indicate that there is 
also a deficit of i/n — > i/n. 

This measurement represents the first attempt at a 
As determination using this ratio. The systematic er- 
rors are quite large, mostly due to large uncertainties 
in the optical model of the mineral oil. MiniBooNE 
maintains a sensitivity to As for proton energies above 
Cherenkov threshold for protons, where the contribution 
from NCE-like background is significant. In order to 
improve the sensitivity to As, future experiments need 
to have good proton/neutron particle identification (pos- 
sibly through neutron capture tagging) and extend the 
cross-section measurement down to the T < 200 MeV re- 
gion, where the contribution from NCE-like background 
becomes negligible and where the extrapolation of the 
axial form factor to Q^ = becomes less model depen- 
dent. 



data sample includes 7,616 events with the following pre- 
dicted channel fractions: 55% vp — >■ i/p, 10% vn — > vn, 
14% NCE-like background, 1% dirt background, and 19% 
other backgrounds. 

The ratio of lyp — > i/p to vN -^ vN events for data 
and MC for different values of As (—0.5, 0, and +0.5) is 
shown in Fig. [TH which illustrates the sensitivity of this 
ratio to As. The error bars for the data histogram are 
the diagonal elements of the full error matrix. 



C. Measurement of As. 

The full error matrix is used for the y^ tests of As to 
determine the best fit and the confidence interval. The y^ 
surface also slightly depends on the value of Ma ■ Assum- 
ing Ma = 1.35 GeV, the fit to the MiniBooNE measured 



V. SUMMARY. 

In summary, MiniBooNE has used a high-statistics 
sample of NCE interactions to measure the NCE [vN -^ 
vN) flux-averaged differential cross section, da/dQ^ on 
CH2. Using MiniBooNE CCQE data, a measurement 
of NCE/CCQE cross-section ratio has also been per- 
formed. Using POT-normalized distributions of the re- 
constructed energy for the NCE sample, x^ tests for sev- 
eral Ma and k values have been performed. The MC 
with higher values of Ma give a better x^ than that with 
Ma = 1.02 GeV. The allowed region for the axial vec- 
tor mass using just MiniBooNE NCE data was obtained. 
Ma = 1.39±0.11 GeV. It is in agreement with the shape 
normalized fits of v^ CCQE scattering on neutrons bound 
in carbon, oxygen and iron as obtained by recent experi- 
ments 0-11, [3 • 1^0^ energies above Cherenkov threshold, 
a sample of NCE proton-enriched events was obtained. 
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which was used for the measurement of the vp — >■ vp to 
vN — ^ vN ratfo, which in turn is sensitive to As. A value 



of As = 0.08 ±0.26 was extracted, in agreement with the 
results from the BNL E734 experiment [13 [. 
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Appendix A: Dirt Background Measurement. 

The fact that dirt and in-tank events have different 
spatial distributions, such as reconstructed radius, R, 
and the reconstructed Z coordinate (which is in the di- 
rection of the beam) can be used to determine the dirt 
energy spectrum. Dirt event vertices are generally re- 
constructed closer to the edge of the detector, than the 
in-tank events. Also, dirt events are mostly reconstructed 
in the upstream part of the detector with Z < 0, whereas 
in-tank events have approximately a uniform distribution 
in this variable. Additionally, we use the fact that dirt 
and in-tank events have a very different energy spectrum 
as seen in Fig[4l 

The dirt energy spectrum is measured by fitting the 
MC in-tank and dirt templates to the data in Z, R, 
and energy distributions for the dirt-enriched event sam- 
ples. The R and Z distributions are obviously correlated. 
However, the samples used for the Z and R fits of the 
dirt have a large fraction of events that are present in 
one sample and not in the other. 

To measure the dirt background in the NCE event sam- 
ple, three additional dirt-enriched samples of events are 
used. For each of the variables (reconstructed Z, R and 
energy) different samples are used, based on cutting on 
the other of these variables; for example, if measuring 
dirt events from the Z distribution, one would have an 
additional cut on R. 

The samples that are selected for these fits are defined 
in Table |TT1 The following precuts are the same for each 
event sample (cuts from (1) through (6) described in Sec- 
tion |TTTB|): 

Precuts = 1 Subevent + VHits < 6 + THits > 24 

+ 4400 ns < Time < 6500 ns + T < 650 McV 
+ \n{CelCp) < 0.42. 

These are the same as the NCE analysis cuts with the 
removal of the radial cut. The dirt fraction in the dirt- 
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enhanced samples is increased significantly, by a factor 
of 2-3 over the unenhanced sample. 



Dirt Rate Measurement from the Reconstructed Z 
Distribution. 



To measure the dirt background from the Z distribu- 
tion, the "Dirt-Z" event sample from Table |ll] is used. 
The shapes (templates) of the Z distribution for the in- 
tank and dirt events in the MC are used to fit the shape 
of the Z distribution for the data. These fits are done 
in bins of reconstructed energy T, so that in the end one 
obtains the measured dirt background energy spectrum. 
From the fits the correction function 



Corr 



Dirt 



(Data^/MC^) 



Dirt 



that should be applied to the MC spectrum of the dirt 
events is extracted, which is the ratio of the measured to 
predicted number of dirt events in the i-th bin. 

An example of the fit is shown in Fig. [T3] for several 
intervals of the reconstructed energy. The agreement be- 
tween the data and MC is much better after the fit. 



Dirt Rate Measurement from the Reconstructed R 
Distribution. 



This procedure is essentially the same as for the dirt 
measurement using the Z distribution, but instead, the 
"Dirt-R" sample from Table [TT] is used and the fitting is 
done for the R variable. Again, from the fits, the cor- 



Dirt 



rection function for the dirt energy spectrum Corr 
extracted. 



Dtrt Rate Measurement from the Reconstructed Energy 
Distribution. 



For this method, two event samples are used, the signal 
(NCE) and the dirt-enriched sample with the "Dirt-E" 
cuts from Table |IT1 For both samples, we look at both 
the MiniBooNE data and the MC prediction for NCE, 
dirt and in-tank backgrounds. 

Assuming that the fractions of signal and dirt events 
in both samples are stable relative to MC variations, one 
can measure the spectrum of dirt events in the NCE sam- 
ple from the data distribution for both of these samples. 
We define the following histograms: 

V — reconstructed energy spectrum for data 

B — reconstructed energy spectrum for MC in-tank 

backgrounds 
S — reconstructed energy spectrum for MC NCE 
D — reconstructed energy spectrum for MC dirt. 



which have upper indices describing the event sample, 
namely 

s — NCE event sample 

d — "Dirt-E" event sample. 

The new spectra for the dirt and signal events can be 
determined from fitting the data in both NCE and "Dirt- 
E" samples. In terms of the definitions that we have 
introduced, the condition that these spectra coincide in 
both event samples can be written as: 




For each reconstructed energy bin i there are 6 unknowns 
on the left hand side and 2 knowns on the right hand side 
(the data in both NCE and "Dirt-E" samples). Assum- 
ing a reliable in-tank background prediction, one can fix 
Bf and Bf. Then, we introduce the fractions of signal 
and dirt events in the two samples: 






and 



Sf 



Because these variables are ratios, they are relatively sta- 
ble to MC variations and independent of the dirt and 
NCE events energy spectra. The functions / and g are 
determined from the MC. 

Herewith, one can express Df (the dirt energy spec- 
trum in the NCE sample) in terms of the above defini- 
tions as: 



D',= 



g, {yt - Bt) -^{yf- Bf) 
9t- fi 



which is the measured spectrum of dirt events in the NCE 
event sample. 

Finally, using all three methods, the combined dirt en- 
ergy spectrum correction function fit is performed. The 
chosen form of the fit function is linear below 300 MeV 
and a constant above, as shown in Fig. [T31 All of these 
measurements agree with each other within 10%, lend- 
ing confidence to our overall ability to constrain the dirt 
background in the analysis. 

The new, measured reconstructed energy spectrum for 
the dirt events in the NCE sample is calculated as a 
bin-by-bin correction of the initial MC dirt energy spec- 
trum multiplied by the measured correction function. 
The measured number of dirt events is determined to 
be ^ 30% lower than the original MC prediction. 



Appendix B: MiniBooNE Neutral— Current Elastic 
Cross-Section Description. 

B.l Phenomenology. 

While MiniBooNE uses relativistic Fermi gas model 
of Smith and Moniz to describe NCE scattering, it is 
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Sample name Purpose of the sample 



Cuts : Precuts+ 



Dirt fraction (%) 



NCE 
Dirt-Z 
Dirt-R 
Dirt-E 



NCE sample (dirt-reduced) Rfiduciai{T) (cut (7) in Section Till B|) 

Fit dirt from Z (dirt-enhanced) 3.8 m < 7? < 5.2 m 

Fit dirt from R (dirt-enhanced) Z < m 

Fit dirt from energy (dirt-enhanced) 3.8 m < i? < 5.2 m and Z < m 



13.4 
27.8 
34.3 
37.6 



TABLE II: Event sample cuts, their respective purposes, and dirt events fractions. The dirt fractions are calculated from the 
initial MC simulation (before dirt fits), NCE is the signal sample (cuts from 1 to 7), whereas the other three samples are the 
dirt-enriched samples for use in the dirt fits. 
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FIG. 13: Fits to the data using MC templates for the in-tank and dirt events in Z for different reconstructed energy bins. The 
"total MC before the fit" is a sum of the in-tank and dirt templates, which are absolutely (POT) normalized. 



illustrative to write down the expression for the NCE 
cross-section in case of a free nucleon target. 

The neutrino-nucleon NCE differential cross-section on 
free nucleons can be written as ill;]: 



where the '-I-' sign corresponds to neutrinos and the '— ' 
sign to anti-neutrinos, W = AE^/M^f — Q'^/M^, and 
A{Q^), B{Q'^) and C{Q'^) are form factors defined as: 



da GJ,Q^ 
dQ^ ~ 2'kEI 



(A(g2) ± B{Q^)W + C{Q^)W^ 



(2) 
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FIG. 14: Combined dirt energy correction fit from Z, R and 
energy distributions. Errors for Z and R are statistical from 
the minimization fitting these distributions. The error for 
the energy distribution is the largest systematic uncertainty 
resulting from the optical model of the mineral oil. 






[(^l)'(l + r)- {{F^f - r{Fif){l - r) + AtF^ F^] 



-/iiFt 



F4 



M 



N 



16Q" 



mr+{Ffr+r{Fir] 



Here Ff , F2 , and F^ are nucleon Dirac, Pauli and axial weak neutral-current form factors, respectively, which in 
general are real dimensionless functions of Q^, and t = Q'^ /AM^. Each of the nucleon for factors is different for proton 
and neutron targets. At low Q^ the C{Q^) term in Eq.Q dominates (see Ref . |3l| ) . Thus, the NCE cross-section has 
a significant contribution from axial vector currents. 

Under the conserved vector current [36|, one can express the weak form factors through their electromagnetic 
equivalents: 



F' 



sm 6w 



F. 



EM,p 



F, 



EM.i 



T3 



1W 



pEALp . pEMj 



-F'' 



1,2, 



(3) 



where 9w is the weak mixing angle, T3 is a factor +1 for 
protons and -1 for neutrons, F^ are the isoscalar form 
factors discussed later, and the Dirac and Pauli electro- 
magnetic form factors are: 



Ge{Q^) + 7?7rGM(g') 



Ff'^Q^) 



Fr'HQ') = 



4M 



N 



Q' 



4M 



N 



GMiQ^)-GE{Q^) 



l2 ' 



4.Mff 



where Ge and Gm are Sachs form factors l37[. Instead 
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of the dipole approximation the BBA-2003 form of the 
Sachs form factors [Sg] is used in this analysis, which 
better describes the electron-proton scattering data. 

The axial weak form factor by definition can be ex- 
pressed through its isovector and isoscalar parts: 



fl 



T3 



Fa 



rFl 



(4) 



The isovector axial form factor can be measured via weak 
charged-current. Usually it is assumed to have a dipole 
form: 



Fa{Q' 



9A 



1 + 



M 



(5) 



where gA — Fa(S^) = 1.2671 [3J] is measured precisely 
from neutron beta decay. 

The isoscalar form factors Ff and F| in Eq.(I3]) are 
usually thought to be due to contributions from strange 
quarks to the electric charge and to the magnetic moment 
of the nucleon, whereas F\ in Eq.(j4]) is their contribu- 
tion to the nucleon spin. F^ value at Q^ = is called 
As. Ff and F2 can be extracted from parity violating 
electron scattering experiments. Recent results from the 
HAPPEX experiment [33! show that these electric and 
magnetic strange form factors are consistent with zero. 
The exact expression for the axial isoscalar form factor 
is unknown, but in analogy to the isovector axial form 
factor it is usually represented in the dipole form with 
the same value of the axial vector mass to minimize the 
number of free parameters in the model [l3| : 



FXiQ^) = 



Ml 



(6) 



B.2 MiniBooNE Neutral-Current Elastic 
Cross-Section Discussion. 

The MiniBooNE NCE scattering sample consists of 
three different processes: scattering on free protons in hy- 
drogen, bound protons in carbon, and bound neutrons in 
carbon. Because several final state nucleons may be pro- 
duced, we define the interaction in carbon using most en- 
ergetic final state nucleon. This means, for example, that 
it is possible for an event to be tagged as a NCE neutron, 
because it has a neutron as the most energetic final state 
nucleon, even though the original neutrino interaction 
was on a proton. According to NUANCE, the probability 
of this misidentification grows almost linearly from about 
8% at Ql^E = O.lGeV^ to 16% at Q|^ = 1.6 GeV^ 

The result shown in Fig. [5] is the flux-averaged NCE 
differential cross-section on CH2 , averaged over these pro- 
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FIG. 15: Efficiency corrections Cvp.H, Cvp,c, and Cvn,c ver- 
sus Qqe for difi'erent NCE processes (as labeled). 



cesses. Herewith, the uN — > vN cross-section is ex- 
pressed as: 



da, 



uN^uN 



1 



dQ' 



■Z Cup^h{Qqe) 
+ T^Cjyp^ciQoE) 



da. 



da, 



wp^up.H 



7 -''P^^WqB^ dQ2 



(7) 



+ -jCunfiiQqE) 



da. 



un^^un^C 



dQ' 



where da. 



Up^!-Vp 



H /dQ^ is the NCE cross-section on free 
protons (per free proton), da^p-^vp,c / dQ'^ is the NCE 
cross-section on bound protons (per bound proton), and 
davn-¥vn,c / dQ"^ is the NCE cross-section on bound neu- 
trons (per bound neutron). The efficiency correction 
functions C^p^H, G^p^Ci and C^n,c result from different 
selection efficiencies for each type of NCE scattering pro- 
cess, and are estimated from the MC as functions of Qg^ 
- as shown in Fig. [121 They are defined as the ratios of 
the efficiency for a particular type of NCE event to the 
average efficiency for all NCE events in bins of Q'qe- 

As one can see from Fig. [151 the efficiency correction 
functions are equal to unity in the region of Qg^ from 
0.4 to 1.2 GeV^ However, at higher Q^^, NCE neu- 
trons have higher efficiency, thus having a higher proba- 
bility than NCE protons to pass the T < 650 MeV cut 
used in the selection. Similarly, at lower Qg^;, NCE neu- 
trons have lower probability than NCE protons to pass a 
THits > 24 cut. 

To calculate a cross-section which is to be compared to 
the MiniBooNE results (as in Fig. [5]), one needs to apply 
these efficiency corrections to each predicted distribution. 
However, note that in the bulk of the measured region, 
0.1 < Q^E < 1.0 GeV", they are all roughly equivalent 
to unity. 

Finally, we present tables with the NCE/CCQE diflFer- 
ential cross-section ratios in Table IIIII and the NCE dif- 
ferential cross-section, NCE-like background and correc- 
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O^iC pV^ ^ \ Distribution 


O-JVCE 
CCCQE 








CNCE-Uke 
(yCCQE-Uke 




0.100-0.150 


(2.019 ±0.261) 


X 


10" 




(1.682 ±0.199) X 10" 




0.150-0.200 


(1.839 ±0.225) 


X 


10" 




(1.582 ±0.181) X 10" 




0.200-0.250 


(1.769 ±0.199) 


X 


10- 




(1.554 ±0.165) X 10" 




0.250-0.300 


(1.696 ±0.173) 


X 


10" 




(1.526 ±0.147) X 10" 




0.300-0.350 


(1.619 ±0.167) 


X 


10" 




(1.502 ±0.145) X 10" 




0.350-0.400 


(1.620 ±0.192) 


X 


10" 




(1.561 ±0.172) X 10" 




0.400-0.450 


(1.594 ±0.231) 


X 


10" 




(1.601 ±0.210) X 10" 




0.450-0.500 


(1.602 ±0.264) 


X 


10" 




(1.680 ± 0.244) X 10" 




0.500-0.600 


(1.584 ±0.313) 


X 


10" 




(1.767 ± 0.295) X 10" 




0.600-0.700 


(1.532 ±0.416) 


X 


10" 




(1.855 ± 0.397) X 10" 




0.700-0.800 


(1.498 ±0.572) 


X 


10" 




(1.966 ±0.554) X 10" 




0.800-1.000 


(1.421 ±0.755) 


X 


10" 




(2.028 ± 0.737) x 10" 




1.000-1.200 


(1.408 ±0.712) 


X 


10" 




(2.102 ±0.720) X 10" 




1.200-1.500 


(1.226 ± 0.643) 


X 


10" 




(1.897 ±0.666) X 10" 




1.500-2.000 


(8.948 ± 5.248) 


X 


10" 


-2 


(1.431 ± 0.736) X 10" 





tion coefficients in Table IIVI There is also an alternative 
method of reporting the results of this analysis, namely 
in terms of the reconstructed nucleon kinetic energy. The 
tables for the latter can be found in 31] and 40] . In or- 
der to make use of these results, one would have to follow 
the instructions described in Appendix B of Ref. j31| . 



TABLE III: NCE/CCQE and NCE-Iike/CCQE-like differen- 
tial cross-section ratios as a function of Qqe = 2m]v X^j Ti. 
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Q2(GeV2)\°^'^"^^*^°'' 


NCE cross-section 




NCE- like background, 


Cvp^H 


Cvp,c 


Cun,C 




cmVGeV^ 






cmVGeV^ 








0.101-0.169 


(3.361 ± 0.360) 


X 


10" 


-39 


4.875 X 


io-« 


0.784 


1.068 


1.014 


0.169-0.236 


(2.951 ± 0.394) 


X 


10" 


-39 


4.623 X 


io-« 


1.206 


1.192 


0.785 


0.236-0.304 


(2.494 ± 0.429) 


X 


10" 


-39 


5.879 X 


io-« 


1.102 


1.119 


0.871 


0.304-0.372 


(2.089 ± 0.340) 


X 


10" 


-39 


8.425 X 


io-« 


1.053 


1.043 


0.949 


0.372-0.439 


(1.744 ± 0.243) 


X 


10- 


-39 


1.235 X 


^0-40 


1.024 


1.009 


0.985 


0.439-0.507 


(1.432 ± 0.246) 


X 


10" 


-39 


1.647 X 


io-*« 


1.016 


1.002 


0.994 


0.507-0.574 


(1.168 ±0.260) 


X 


10" 


-39 


1.964 X 


io-« 


0.986 


0.994 


1.009 


0.574-0.642 


(9.435 ± 2.400) 


X 


10" 


-40 


2.155 X 


io-« 


1.007 


0.989 


1.008 


0.642-0.709 


(7.534 ± 2.205) 


X 


10" 


-40 


2.222 X 


io-« 


1.023 


1.002 


0.992 


0.709-0.777 


(6.015 ±2.194) 


X 


10" 


-40 


2.215 X 


io-« 


1.023 


0.995 


0.998 


0.777-0.844 


(4.832 ± 2.320) 


X 


10" 


-40 


2.075 X 


io-« 


0.994 


1.003 


0.999 


0.844-0.912 


(3.854 ±2.331) 


X 


10" 


-40 


1.890 X 


10-40 


0.993 


1.009 


0.994 


0.912-0.980 


(3.209 ± 2.330) 


X 


10" 


-40 


1.756 X 


io-« 


0.999 


1.004 


0.997 


0.980-1.047 


(2.649 ±2.117) 


X 


10" 


-40 


1.521 X 


10-40 


1.074 


0.980 


0.997 


1.047-1.115 


(2.226 ±1.818) 


X 


10" 


-40 


1.265 X 


10-40 


1.056 


1.015 


0.973 


1.115-1.182 


(1.935 ±1.295) 


X 


10" 


-40 


1.081 X 


10-40 


0.913 


1.034 


0.994 


1.182-1.250 


(1.598 ±0.939) 


X 


10" 


-40 


9.324 X 


10-41 


0.731 


0.967 


1.099 


1.250-1.317 


(1.329 ± 0.769) 


X 


10" 


-40 


8.079 X 


10-41 


0.716 


0.867 


1.187 


1.317-1.385 


(1.111 ±0.689) 


X 


10" 


-40 


6.983 X 


10-41 


0.567 


0.836 


1.251 


1.385-1.452 


(9.259 ± 6.254) 


X 


10" 


-41 


5.958 X 


10-41 


0.504 


0.808 


1.312 


1.452-1.520 


(7.975 ± 5.373) 


X 


10" 


-41 


5.229 X 


io-« 


0.419 


0.832 


1.286 


1.520-1.588 


(6.618 ±4.524) 


X 


10" 


-41 


4.342 X 


io-« 


0.378 


0.746 


1.363 


1.588-1.655 


(5.799 ±3.921) 


X 


10" 


-41 


3.839 X 


10-41 


0.410 


0.818 


1.309 



TABLE IV: MiniBooNE measured NCE differential cross-section, predicted NCE-like background, and predicted correction 
coefficients for the three different NCE scattering contributions as a function of Qqe ~ 2mN X^j Ti. 



